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We use density functional theory (DFT) to investigate the electronic structure and chemical 
properties of gold nanoparticles. Different structural families of clusters are compared. For up 
to 60 atoms we optimize structures using DFT-based simulated annealing. Cluster geometries are 
found to distort considerably, creating large band gaps at the Fermi level. For up to 200 atoms we 
consider structures generated with a simple EMT potential and clusters based on cuboctahedra and 
icosahedra. All types of cluster geometry exhibit jellium-like electronic shell structure. We calculate 
adsorption energies of several atoms on the cuboctahedral clusters. Adsorption energies are found 
to vary abruptly at magic numbers. Using a Newns- Anderson model we find that the effect of magic 
numbers on adsorption energy can be understood from the location of adsorbate-induced states with 
respect to the cluster Fermi level. 

PACS numbers: 36.40.Cg, 36.40.Mr, 36.40.Jn, 34.35.+a 



I. INTRODUCTION 

A major theme in advanced materials design today is 
the possibility to modify and change materials properties 
through structuring at the nanoscale. The applications 
can be as diverse as optimizing the size of metal nanopar- 
ticles to catalyze certain chemical reactions^ or the struc- 
turing of surfaces and interfaces for optimal light absorp- 
tion in photovoltaic devices^. In very broad terms the 
interesting possibilities arise when the structures reach 
a scale comparable to the wavelengths of the relevant 
quantum particles (electrons, plasmons or photons). 

In this work we investigate theoretically the proper- 
ties of freestanding metal nanoparticles made of gold in 
particular. The purpose is to improve our understanding 
of the relationship between cluster size and a range of 
electronic and chemical properties. Different aspects of 
this has been investigated in numerous studies. See for 
example the review by Baletto and Ferrando.^ What is 
special here is that we investigate the cluster properties 
over an — for electronic structure calculations — unusually 
large size range and for many different cluster structures. 
The hope is thereby to get a more complete picture of the 
general trends in the cluster behavior. 

For transition metals with partially filled d-bands, co- 
hesive energies will be dominated by the effect of the d- 
states^. Because of the short range of the d-states, their 
contribution to the cluster energy is determined mostly 
by the local arrangement of neighboring atoms. Facet 
types and local atomic packing can therefore be expected 
to be particularly important factors in the structures of 
transition metal clusters with partially filled d-bands. 

The effect of the partially occupied d-band disappears 
for noble metals and alkali metals. Instead the long-range 
s-electrons, which hybridize in a more complex manner, 
yield the primary contribution to the cluster energy. The 
optimal structure will not be determined by optimizing 



the local structure around each atom, but rather by opti- 
mizing the global geometric structure to obtain the most 
desirable electronic structure of the delocalized electron 
cloud. The result is a much more complicated interplay 
between electronic and geometric structure. 

Small free-standing gold clusters have been theoreti- 
cally shown to possess very diverse ground-state geome- 
tries depending on cluster size. Examples are planar, 
cage-like and tube-like structures'^. 

The s-electron hybridization can be interpreted in 
terms of a jellium model which regards the whole clus- 
ter as a spherical superatom. The s electrons organize 
into global shells, resulting in electronic "magic num- 
bers" when shells are filled. Magic numbers at 2, 8, 18, 
20, 34, 40, 58, . . . , have been observed as particularly 
stable alkali metal clusters^ii^ with large band gaps in 
agreement with theory. For alkali metals, magic num- 
bers attributed to both electronic and geometric shell 
structures have been observed for clusters with thousands 
of atoms, with geometric shells dominating beyond 2000 
atoms^^. Larger Au clusters are believed to form icosa- 
hedra, decahedra or truncated octahedra depending on 
size and temperature^^"— . 

In this work we consider several series of clusters based 
on different generation procedures and structural motifs. 
We calculate structures of smaller clusters using simu- 
lated annealing with density functional theory (DFT) 
and for larger clusters using effective medium theor}^!^^ 
(EMT). Using DFT we compare the energy and elec- 
tronic structure of optimized clusters with the commonly 
considered regular icosahedral and cuboctahedral struc- 
tures. For the cuboctahedra we identify trends in reactiv- 
ity by considering adsorption of different atoms. The ge- 
ometric similarity of clusters based on cuboctahedra and 
icosahedra allows us to isolate and study size-dependent 
effects on chemistry. The price of this simplification is 
that individual calculations do not represent globally op- 
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timal structures. Hence we focus on trends that are gen- 
eral enough to be significant outside the model systems. 

II. COMPUTATIONAL METHODS 

All electronic structure calculations are performed 
with the real-space DFT code GPAW^^^^^ using the 
RPBE^ functional for exchange and correlation. GPAW 
uses the projector augmented wave (PAW) method^^, 
and offers an accurate real-space representation of the 
Kohn-Sham orbitals along with an efficient basis set of 
localized atomic orbitals^-. The calculations presented 
here are performed with the atomic orbitals using a 
double-zeta polarized (DZP) basis set. All calculations 
on clusters are spin-paired and use the scalar-relativistic 
atomic PAW setups and basis sets supplied with GPAW. 
The Au setup contains 11 valence electrons. 

In our calculations the cluster is centered in a non- 
periodic orthorhombic cell with 5.0 A vacuum along each 
axis. We use a grid-spacing of 0.2 A. 

We do not apply any basis set superposition error cor- 
rection, so the values of adsorption energies are not neces- 
sarily accurate. However in comparing the bonding of an 
adsorbate to clusters of different sizes, the local structure 
around the adsorbate highly similar for all clusters, and 
the basis set error is consequently roughly the same for 
all clusters. Therefore variations in adsorption energies 
are subject to a much smaller error. 

Pulay density mixing^^ is used to speed up convergence 
of the self-consistency loop. Electron occupations are 
smeared by a small Fermi temperature of 0.01 eV, which 
helps speed up convergence. 

Structure optimizations are performed using the BEGS 
algorithm as implemented in the Atomic Simulation 
Environment^^ (ASE), and terminate when no force on 
any atom is larger than 0.075 eV/A. 

III. CLUSTER GEOMETRY 

Systematic calculation of lowest-energy structures 
from first principles is computationally very expensive. 
Previous studies of structures and properties of Au clus- 
ters have therefore usually been limited to a few dozens 
of atoms Here we focus mainly on larger clusters 
which are quite challenging to systematically optimize 
and characterize, but which are clearly of interest both 
from a conceptual point of view and in applications like 
catalysis. The transition from smaller clusters over larger 
clusters to bulk-like behavior has been studied recently 
by Kleis et al?^ and the results presented here can be 
seen as a supplement and expansion of this study. 

In this work we compare clusters generated by several 
different procedures. For the smallest clusters we perform 
simulated annealings using DFT to obtain realistic struc- 
tures. This is clearly the most realistic and theoretically 
satisfactory method since the same energy landscape is 



used to define the cluster geometry or shape as is used 
to subsequently study the bonding of the cluster and the 
chemical properties. However, for reasons of computer 
time this approach cannot be generally applied for larger 
clusters. 

Larger clusters are studied using DFT, but with the 
structures being determined by simulated annealing with 
a classical EMT potential. As this EMT potential does 
not incorporate explicit electronic structure, the struc- 
tures generated by this method will have no information 
about potential electronic shell effects but only of atomic 
shell effects related to atomic packing of the clusters. 

We finally construct clusters based on prescribed 
cuboctahedral and icosahedral shapes. The simplicity of 
the fcc-based cuboctahedral structures allows us to study 
adsorption of atoms in a way which preserves the local 
geometry around the adsorbate for different cluster sizes. 
This allows us to separate the effect of local geometry 
from that of the electronic structure of the cluster, which 
would not generally be possible if the cluster were based 
on a global minimum search. The comparison of distinct 
types of structures will help determine how properties of 
clusters depend on structure versus size. 



A. Simulated annealing with DFT 

For the smallest clusters (A/'=6-60) we calculate real- 
istic geometries using DFT with coarse parameters. 

For each size of cluster we perform a rough simulated 
annealing based on molecular dynamics (MD) to find the 
optimal structure. We use a Langevin thermostat to reg- 
ulate the temperature from 750 K to 300 K. For a clus- 
ter of size TV we lower the temperature by 1 K for every 
5 + N/2 timesteps of length 24 fs. The timestep is too 
large to have accurate energy conservation during the op- 
timization. This can cause unrealistic behavior when the 
atoms move quickly, but is not likely to affect the results 
of an annealing where the result is mostly determined at 
lower temperatures. The optimization is performed with 
a very coarse grid spacing of 0.24 A. 

At the end of the MD simulation we perform a struc- 
ture optimization with normal DFT parameters using the 
BEGS algorithm such that the structure is guaranteed to 
be a local minimum. 

These optimizations produce planar and tetrahe- 
dral structures in qualitative agreement with previous 
finding9^i2^"— , while larger structures tend to be irregu- 
lar but with some well- formed facets. Due to the short 
annealing times, the larger structures are unlikely to be 
global optima. Figure [1] shows the 20-atom tetrahedron 
and the 58-atom cluster obtained with this method. No- 
tice on the 58-atom cluster the imperfect five-fold sym- 
metry center reminiscent of those found on icosahedral 
clusters. 



FIG. 1: (Color online) 20 (left) and 58- atom (right) clusters 
obtained by simulated annealing with DFT. 




FIG. 2: (Color online) Cuboctahedra (top) and icosahedra 
(bottom) with 55, 147 and 309 atoms. For the cuboctahedra, 
an O atom is shown at the fee site closest to the center of an 
111 facet. 



B. Simulated annealing with EMT 

For a larger range of clusters (A/'=6-200) we calcu- 
late structures using a simple EMT potential^^^^^ imple- 
mented in AS AP^^ . This potential is designed to provide 
reasonable descriptions of elastic and cohesive properties. 
It is a classical potential and as such contains no explicit 
description of electronic behavior. 

For each size of cluster we perform a simulated anneal- 
ing wherein the temperature varies from above the melt- 
ing point (1337 K for Au) to 200 K with 600 MD steps of 
6.0 fs for each IK decrease in temperature. The result- 
ing structures frequently have 5-fold symmetry centers 
surrounded by 111 facets, resembling partially formed 
icosahedra or decahedra. Again, at the end of the simu- 
lated annealing we perform a BFGS structure optimiza- 
tion with the usual DFT parameters. 



C. Atomic shell structures 

Finally we consider the cuboctahedral and icosahedral 
series of structures. Each cluster can be constructed geo- 
metrically from the previous one by adding one complete 
shell of atoms. The cuboctahedra and icosahedra have 




FIG. 3: (Color online) Generation of regular clusters with 
different numbers of atoms. The white atoms belong to the 
55-atom cuboctahedron, while the grey atoms are stripped off 
one by one as marked with a cross. Removable atoms with 
lower coordination numbers have darker shades of grey, and 
at each step one of the lowest- coordinated atoms is removed 
at random. An oxygen atom is shown at the adsorption site. 



closed atomic shells at the same numbers. The first few 
geometric shell closings are A^=13, 55, 147, 309 and 561. 

We would like to study the chemical properties of clus- 
ters by calculating adsorption energies of atoms on clus- 
ters of different sizes. A systematic comparison can be 
made if we ensure that the local geometry around the 
adsorbate remains identical for all clusters independent 
of size. 

For the cuboctahedra we generate clusters with dif- 
ferent numbers of atoms by stripping off atoms one by 
one from one cuboctahedron until only the next smaller 
cuboctahedron remains. For each step, the next atom 
to be removed is chosen at random amongst those that 
have the lowest coordination numbers and are not part 
of the smaller cuboctahedron. This procedure is shown 
on Figure [3l 

As mentioned, to obtain adsorption energies that can 
be meaningfully compared across the different clusters, 
we must avoid changing the immediate environment of 
the adsorbate when removing atoms. For this reason we 
do not simply remove the outermost shell. Instead we 
choose an adsorption site on a particular facet, then re- 
move atoms as necessary on the other sides of the cluster 
such that the local facet is changed only minimally. 

This method preferentially strips off corner atoms and 
atoms on the most open facet, opening a new facet only 
when necessary. This avoids very unphysical geometries. 
We run the calculations multiple times using a pseudo- 
random number generator with different random seeds. 
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FIG. 4: (Color online) Energy per atom for different Au clus- 
ter types as a function of cluster size. The energy reference 
corresponds to bulk Au. 



yielding a small ensemble that shows the dependence of 
cluster properties on the randomization. 

A similar procedure can be applied to the icosahedra. 
However in the icosahedra, the distance between atoms 
in successive layers is different from the distance between 
atoms within the same layer, which means the local ge- 
ometry around an adsorbate cannot always be preserved 
as for the cuboctahedra. 



IV. INTRINSIC PROPERTIES OF CLUSTERS 

In the following we study energies and electronic prop- 
erties of the different types of clusters. Note that all 
clusters have been relaxed using DFT such that all struc- 
tures are local minima corresponding to the same force 
method, and have directly comparable total energies. 



A. Cluster structure and stability 

Figure S] compares the energy per atom for Au clus- 
ters of the different types. EMT-optimized and reg- 
ular clusters have been generated multiple times from 
different pseudorandom number sequences, yielding four 
datapoints for each cluster size. The DFT-optimized 
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FIG. 5: (Color online) Energy for different Au cluster types 
as a function of cluster size N. A smooth function of N cf. 
Eq. [T]has been subtracted from all datapoints for legibility. 



structures generally have energies lower than or equal 
to the other methods followed by EMT. Among the reg- 
ular structures, icosahedra usually have lower energies 
than cuboctahedra. Even where the regular clusters have 
closed geometric shells (A/'=55 and 147), they are less sta- 
ble than the structures obtained by simulated annealing 
with EMT. 

Prominent kinks in the energy are visible around 
A/'=34, 58 and 92 atoms. These are "magic numbers" 
corresponding to major electronic shell closings. They 
are well known in the jellium cluster models of simple 
metals^^^^, and have also been observed in mass spectra 
of noble metal clusters^^. The kinks in energy due to 
electronic shell structure are robust enough to be visible 
for all types of clusters considered. Figure [5] provides a 
closer view of the energies of smaller clusters. To improve 
legibility, a smooth function of N of the form 

^fit(^') = ao + aiN^^^ + asA/-'/^ + asN, (1) 

is subtracted from all energies. The coefficients 
(ao, ai, a2, as) are obtained by fitting the energies of the 
DFT-optimized clusters. For the other kinds of clusters, 
only the lowest-energy datapoint found among four at- 
tempts is shown for each N. 

The DFT-optimized clusters up to 13 atoms are planar 
except for A^=10 and 11. The predicted transition be- 
tween planar and three-dimensional structures depends 
strongly on the approximation of exchange and correla- 
tion, and has been studied more systematically by several 
authors^^^^^. 

The particularly visible feature at A^=20 is the well- 
known tetrahedron^^. Aside from this, particularly sta- 
ble clusters are Ar=34, 40, 48 and 56/58. The EMT- 
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FIG. 6: (Color online) Square roots of ratios /1//3 and I2/I3 
of the three principal moments of inertia Ii < I2 ^ H as a 
function of cluster size, showing deformations of the struc- 
tures. 



based and regular structures tend to obtain comparable 
energies only around the major magic numbers A/'=34 
and 58. In between the magic numbers, the EMT-based 
structures have higher energies than the DFT-optimized 
ones by typically 1-2 eV. 

The main structural difference between the DFT- 
optimized clusters and other types is that the DFT- 
optimized clusters systematically deviate from spheri- 
cal shapes when doing so is favorable to the electronic 
structure. A rough measure of how spherical a cluster 
is can be obtained by considering the moments of in- 
ertia. For each cluster the three principal moments of 
inertia /i < h ^ h ar e calc ulated. Figure [6] shows 
the ratio y/i/Ta and ^JhjH^ i-e. the square root of 
the inverse ratio between the largest principal moment 
and each of the smaller ones, as a function of cluster 
size. The most symmetric clusters are found around the 
magic numbers A/'=20, 34, 40 and 58, while intermedi- 
ate clusters are deformed considerably. A similar varia- 
tion has been predicted for Cu clusters^ii^^ and in sev- 
eral deformable-background jellium models^SdU. Clus- 
ters just above magic numbers are from jellium models 
expected to be prolate while clusters below magic num- 
bers are expected to be oblate. Such a trend is not clearly 
visible from our results. This is most likely due to the 
roughness of the optimization procedure combined with 
the presence of a physical atomic lattice, modifying the 
simple model picture. 

Figure [7] compares the stabilities of Cu, Ag, and Au 
(top). For comparison two other transition metals Pd 
and Pt are also shown (bottom). For each species, the 
energy is calculated using simulated annealing with EMT 
followed by a local geometry optimization with DFT for 
A^=6-200. A smooth function is then subtracted by fit- 
ting the energies for each element according to Eq. [TJ 
such that the figure shows the deviation from a smooth 
trend line. 
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FIG. 7: (Color online) Energy minus fitted trend line for noble 
metal clusters (top) and d-band metals (bottom). Electronic 
magic numbers at A/'=18, 34, 58, 92 and 138 are indicated. 



As for Au, the other noble metal clusters are par- 
ticularly stable close to the magic numbers 18, 34, 58 
and 92, and to a lesser extent at 138. Deviations from 
the trend line oscillate with a peak-to-peak variation of 
around 3eV. Beyond 138, the periodic trend is gradually 
obscured by fluctuations which are probably caused by 
imperfections in the optimization procedure. 

The three noble metals exhibit roughly identical be- 
haviour except in the region N < 18, where Au differs 
noticeably from Cu and Ag. In this case the Au clusters 
deform considerably from the geometries found by EMT, 
tending towards flat structures. Cu and Ag clusters re- 
main round. The tendency of small Au clusters to form 
planar structures has been well documented and has been 
attributed to relativistic effects causing a contraction of 
the s-orbitals compared to the d-orbitals. — — The differ- 
ences in behavior between noble metals here appear to 
be caused exclusively by the relativistic behavior of Au, 
and not e.g. the location of the d-band in relation to the 
Fermi level, which would set Cu and Ag apart. 

For Pd and Pt, no magic numbers are observed. De- 
viations from the trend line instead appear to depend on 
how well-formed the clusters are, i.e. the type and regu- 
larity of their facets as seen from visual inspection of the 
cluster structures. Thus the stability of non-noble transi- 
tion metal clusters is determined mostly by local atomic 
arrangement, corresponding to interactions between the 
short-ranged d-electrons. 
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FIG. 8: (Color online) DOS in arbitrary units of EMT- 
optimized Au clusters (top) and Pt clusters (bottom) as a 
function of cluster size and energy. The line indicates the 
Fermi level. Values larger than 1.0 are truncated to 1.0. 



B. Electronic structure 

Figure [9] compares the density of states (DOS) per 
atom of Au (top) and Pt (bottom) clusters optimized 
with EMT as a 2-variable function of cluster size and 
energy. For each cluster, the DOS is approximated as 
a sum of Gaussians of width 0.07 eV centered on each 
energy eigenvalue. 

For both Au and Pt the d-states very quickly form 
the usual continuous, narrow band which beyond A^=20 
changes only very little. The s-states split up into multi- 
ple electronic shells which are separated by gaps as in the 
jellium shell model. As N increases the shells gradually 
broaden to form a continuous band. Oscillations in the 
DOS originating from the shell structure are still clearly 
visible even for the largest clusters. 

For Au, where the Fermi level is located well above 
the d-band, the electronic shells due to the s-electrons 
are filled one by one as cluster size increases. When one 
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FIG. 9: (Color onUne) DOS in arbitrary units of DFT- 
optimized clusters. The line indicates the Fermi level. DOS 
values larger than 1.0 are truncated to 1.0. 



shell is full, the Fermi level jumps to the next higher shell, 
causing the abrupt shifts in Fermi level and large band 
gaps. 

The Fermi level for Pt is lodged at the top of the d- 
band where the DOS is extremely high. Therefore no 
gaps or jumps in the Fermi level are possible, and the 
cluster will not exhibit any electronic magic numbers 
even though the s-electrons form shells exactly like Au. 

Comparing to the DOS of the DFT-optimized clusters 
on Fig. [9l the DFT-based optimization consistently cre- 
ates very large gaps around the Fermi level for all small 
clusters. For clusters with an odd number of electrons 
(where the gap is zero because of spin-degeneracy), a 
single singly-occupied state is located at the middle of a 
symmetric gap. Similar behavior has been reported for 
Cu cluster^. A significant difference compared to the 
roughly spherical clusters obtained with EMT is that the 
shell structure cannot easily be seen as distinct shells that 
move down in energy as the cluster size increases. Rather 
there is an accumulation of states some way above as well 
as below the Fermi level. Only close to the shell closings 
does the DOS resemble that of the EMT-optimized clus- 
ters. A consequence of this is that the abrupt jumps in 
Fermi level cf seen for EMT-optimized clusters are less 
visible for the DFT-optimized ones. However a signifi- 
cant change in ionization potential / and electron affinity 
A still accompanies a magic number. Figure [TOk shows 
the difference —{I-\-A)/2 ^ ep- The value increases 
sharply at each magic number. 

It is easy to understand that gaps at the Fermi level 
are associated with an increase in stability. The gap is 
created so that occupied states are pushed down in en- 
ergy while pushing unoccupied ones up, resulting in a 
decrease of band structure energy. 

For larger clusters that are not close to magic num- 
bers (e.g. N ^ 45) the gap becomes small, but a signif- 
icant depletion of states around the Fermi level persists 
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FIG. 10: (Color online) Fermi level (top) and chemical hard- 
ness (bottom) calculated from ionization potential / and elec- 
tron affinity A for DFT-optimized Au clusters as a function 
of cluster size. 



FIG. 11: (Color online) DOS near Fermilevel for different 
types of clusters, (a) Clusters optimized using the EMT po- 
tential, (b) Clusters based on icosahedra. (c) Clusters based 
on cuboctahedra. Units are arbitrary. 



(a similar depletion of states close to the Fermi level is 
also seen for the EMT-optimized clusters e.g. for N ^ 45 
and 70. This is a product of the local structure optimiza- 
tion with DFT after the EMT-annealing). The combined 
structural and electronic trends of the DFT-optimized 
clusters thus point to a picture where clusters far from 
magic numbers will deform significantly, maximizing the 
gap at the Fermi level. In a sense this deformation creates 
a new magic number for every size of cluster provided 
the clusters are small enough. As long as such a gap 
remains, strong even-odd oscillations of the electronic 
properties will persist due to the singly-filled state in un- 
even clusters. Figure [TOb shows the band gap calculated 
as (/ — A)/2 as a function of cluster size. Even and odd 
clusters are plotted as separate lines. The structure opti- 
mization tends to obtain larger gaps close to the spherical 
shell closings, and so the even-odd alternations are larger 
close to these. The even-odd alternations become small 
compared to the 0.1 eV smearing for clusters larger than 
~ 40 atoms except at the electronic shell closings. 

Figure [TTI compares the calculated DOS near the Fermi 
level for EMT-optimized, icosahedral and cuboctahedral 
structures. The structures yield remarkably similar elec- 
tronic shells separated by gaps. 

Similarities between the electronic shell structures of 
spherical and faceted structures with hundreds of elec- 
trons have previously been found in the context of well 
potentials. Our results show that the inclusion of a 
d-band, as well as the inclusion of an atomic lattice with 
various irregularies as per the different cluster genera- 



tion procedures used here have limited effect on the shell 
structure in this size range. 

The highly visible change for icosahedra close to 
A^=130 happens when a sufficiently large number of 
atoms have been removed from the same side of the clus- 
ter, causing a substantial collective movement of the sur- 
rounding atoms (this is therefore just an artifact of the 
cluster generation method). 



V. CHEMICAL PROPERTIES OF CLUSTERS 

In this section we consider the chemisorption of various 
atoms on cuboctahedral clusters. 

Adsorption energies are calculated as follows. A series 
of cuboctahedral clusters is generated by the procedure 
described in Section IIII C( so as to preserve the local ge- 
ometry around the desired adsorption site. A structure 
optimization is then performed on the entire cluster plus 
adsorbate, yielding a total energy of the combined sys- 
tem. From this energy we subtract the energies of the 
isolated adsorbate and the isolated cluster. Ideally the 
energy of the isolated cluster should be calculated by 
removing the adsorbate, then re-relaxing the structure. 
However this is likely to cause atoms to move signifi- 
cantly for at least some of the smaller clusters, obscuring 
the general trends that we are trying to examine with- 
out providing any insight. For this reason we instead 
neglect to re-relax the clusters after removing the adsor- 
bate. Calculated binding energies therefore tend to be 
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FIG. 12: (Color online) Adsorption energies as a function of 
Au cluster size, (a) O on the fee site closest to the center of 
an (111) facet and the hep site at the corner of an (111) facet, 
(b) F on the central (111) fee site, (c) H and Li 



overestimated. 

Figure [T2k shows the adsorption energy of oxygen on 
cuboctahedral clusters as a function of cluster size for 
two different adsorption sites. One is the fee site as close 
as possible to the center of an (111) facet, which locally 
resembles a (111) surface. The other is the hep site clos- 
est to the corner of a (111) facet, where O frequently 
binds more strongly. For each site there are four series 
of datapoints corresponding to different random seeds in 
the cluster generation procedure. 

Adsorption energies on both sites are related to the dis- 
tribution of the electronic shell closings 7V=34, 58, 92 and 
138 where binding energies are particularly low. Near the 
geometric shell closings 55 and 147, where the clusters 
are regular and closer to being spherical, this behaviour 
is most pronounced. The change near 92 is less abrupt, 
and we attribute this to the less symmetric structures 
generated far from geometric shell closings (we believe 
that 92 would stand out more clearly for clusters with 
more realistic geometry). 

Consider the behavior at the magic number 138. The 



binding gradually weakens until 138, after which it 
abruptly changes from very weak to very strong. The 
same effect is seen to a smaller extent at 58 (the preced- 
ing weakening of binding energy is in this case not grad- 
ual, but coincides with the completion of a local facet as 
discussed below). 

Clusters slightly larger than a magic number will have 
one or more loosely bound electrons which can easily be 
donated to O. Clusters can in this sense be characterized 
as alkali-like, noble or halogen-like depending on their 
location relative to magic numbers. What the compar- 
ison between the two adsorption sites shows is that the 
main variation of adsorption energy due to the electronic 
shell structure is not strongly affected by local geometry 
around the adsorbate. While there are intriguing dif- 
ferences between the binding on the central site and the 
corner site, such details are probably too specific to make 
predictions about more realistic geometries. 

The alkali-like or halogen-like behavior of clusters near 
magic numbers is demonstrated on Figures [T2b and [T2b 
which show the adsorption energy of F, H and Li on the 
central (111) fee site as a function of cluster size. 

The variation of F adsorption energy (Figure [T2b) 
around magic numbers is qualitatively similar to that of 
O. Since F is more electronegative, the increase in bind- 
ing energy past a magic number is more abrupt and is 
clearly visible for all the magic numbers 34, 58, 92 and 
138. However F can accept only one electron, and so 
the total magnitude of the increase in energy near A^=58 
and 138 (0.5 eV) is smaller than in the case of O (up to 
1.0 eV). The variation of F binding energy at magic num- 
bers is roughly equal to the increase in Fermi energy of 
the cluster (0.6 eV). 

The electropositive Li shows the opposite trend: a 
steep decrease in binding energy follows a magic number. 
Again, the change in binding energy is close to 0.5 eV 
corresponding to the sharing of one electron. Hydrogen 
somewhat surprisingly follows the same trend as Li even 
though the H Is-state is approximately as low- lying as 
the O 2p-states. We shall analyze this further in the 
next section. 

Apart from the variation due to magic numbers, bind- 
ing energies of all species tend to be stronger for the 
smallest clusters (A^ < 50). The very steep change in 
binding energy around A^=50 which is seen for all ad- 
sorbates can be attributed to geometric changes of the 
local facet. The triangular 6-atom (111) facet of the 55- 
atom cuboctahedron appears to be generally unreactive, 
as all four adsorbates bind weakly to it, including O on 
both the central site and the corner site. Several previ- 
ous works have noted that a lower overall coordination 
of nearby Au atoms has been found to increase binding 
strengths The effects due to local geometry are 
however comparatively small for Au clusters larger than 
55 atoms. The global electronic shell structure is respon- 
sible for most of the variation in adsorption energy, as 
evidenced by the consistent oscillating trend. 

For comparison, adsorption energies of O on Pt clus- 
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FIG. 13: (Color online) Adsorption energy of O on Pt as a 
function of cluster size. 
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FIG. 14: (Color online) Projected density of states on O, F, 
H and Li adsorbed on Au clusters as a function of cluster size 
and energy calculated using the atomic basis set. The white 
line indicates the Fermi level. Units are arbitrary. 

ters are shown on Figure [131 The clusters have the same 
initial structures as the Au clusters in Figure [T2k . In- 
stead of a smooth oscillating trend, the binding energy 
can vary significantly when nearby facets are modified, 
even though the modification takes place several sites 
away. This causes the adsorption energy to depend much 
more sensitively on the geometry of the nearby facets. 
Apart from the strongly geometry-dependent variations, 
an overall decrease in O binding energy with size is also 
observed, which resembles that of Au. We note that the 
changes in adsorption energy on Au clusters occur even 
when atoms are added on the directly opposite side of 
the cluster as seen from the adsorption site. 

Let us return to the Au clusters. The projected density 
of states (PDOS) on the adsorbate reveals useful infor- 
mation on how the atomic states hybridize with the clus- 
ters. The PDOS as a function of energy e on atom A is 
calculated from the atomic expansion of the Kohn-Sham 



eigenstates \ipn) and eigenvalues as 

PAie) = J2{i^n\PA\^n)S{e-en), (2) 

n 

where 

Pa= Yl \^a)iS^%b{^b\ (3) 

aeA,beA 

is a projection operator onto the subspace spanned by the 
basis functions \^a) on atom A, and {SA)ab = {^a\^b) is 
the overlap matrix within that subspace. 

Figure [TH shows the PDOS on the O, F, H and Li atoms 
adsorbed on Au cuboctahedra of varying size. For O the 
p-states split into bonding and antibonding states at the 
top and bottom of the d-band (see also Ref. 47 on hy- 
bridization on Au surfaces). Both bonding and antibond- 
ing states are occupied, and they remain qualitatively 
similar for all sizes of clusters although some variation is 
seen around the magic numbers. The weights of bond- 
ing and antibonding states can change drastically when 
the local geometry is modified, but for clusters larger 
than 55 atoms these changes are not reflected strongly 
in the adsorption energy and do not explain the trends 
(for example the most visible change, around N = 105, 
occurs when the last two second-nearest neighbor atoms 
are added. This causes the adsorption energy to change 
by 0.15 eV). The behavior of F resembles that of O ex- 
cept the coupling is weaker. For H a very low- lying state 
is created at the bottom of the s-band. The surprising 
formation of such a low-lying state at the bottom of the 
band has been described previously for H adsorption on 
Mg surfaces. Finally Li has a high-lying state which is 
above the Fermi level. 



VI. NEWNS-ANDERSON MODEL 

The peculiar properties of metal nanoparticles are 
sometimes attributed to the discrete spectrum causing 
the nanoparticle to behave like a molecule rather than 
a bulk material characterized by a continuous spectrum. 
However the DOS quickly {N > a few dozen atoms) be- 
comes effectively continuous on any reasonable energy 
scale (~ 0.1 eV). The primary feature distinguishing clus- 
ters with a few dozens to a few hundred atoms compared 
to bulk is not whether the DOS is discrete or continu- 
ous, but rather the fact that the approximately continu- 
ous DOS remains grouped into shells separated by gaps. 
The size-dependent variation in this effectively contin- 
uous DOS, and in particular the distribution of magic 
numbers, are the significant factors that make clusters 
with many hundreds of atoms still differ from bulk Au. 
Since the DOS of an Au cluster larger than a few dozen 
atoms can effectively be regarded as continuous, we will 
in the following apply the Newns- Anderson model^^ to 
understand chemisorption on Au clusters. The Newns- 
Anderson model considers a single state \a) on an atom 
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which hybridizes with a continuum of states {\k)} of the 
metal surface described by a Hamiltonian of the form 



H = Hn 



V, 



(4) 



where Hq is the Hamiltonian of the uncoupled cluster 
plus adsorbate, and V describes the coupling. In the 
basis consisting of metal eigenstates {\k)} plus adsorbate 
eigenstate |a), the uncoupled Hamiltonian Hq is diagonal, 
and the Newns- Anderson Hamiltonian can be written in 
block form as 



H 























• • • Vak • • • 


ea - 



(5) 



Here and Ck are the uncoupled energy eigenvalues on 
the atom and in the metal while Vak — (a|V|/c) are cou- 
pling matrix elements. 

Our idea now is to perform a DFT calculation to ob- 
tain a Hamiltonian matrix of the composite system con- 
sisting of both cluster and adsorbate, then transform it to 
Newns- Anderson form cf. Eq. (|5]). In the basis of atomic 
orbit als used by GPAW, the Kohn-Sham equations are 
solved as a generalized eigenvalue problem^ ^ 



(6) 



The overlap matrix S^jy is present because the basis is 
non-orthogonal. The Hamiltonian can be regarded as a 
blocked matrix 



H 



DFT 



j^AM j^A 



(7) 



with one block corresponding to the metal or clus- 
ter, one block H"^ corresponding to the atom and the 
off-diagonal blocks corresponding to the in- 

teraction. We diagonalize each of the metallic and atomic 
submatrices according to 



EttM _ \^ qM m ^ 
^k'v^vk — / , "^k'v^vk^k-, 

V V 



(8) 
(9) 



to obtain values for the energies Ck and of the uncou- 
pled systems. Clearly these are approximate, as the real 
energy values could be evaluated selfconsistently on each 
of the uncoupled systems using a separate DFT calcula- 
tion. However the Hamiltonian and eigenvalues from one 
selfconsistent calculation cannot be expected to be "com- 
patible" with those from a different selfconsistent calcu- 
lation. Indeed it is known from the force theorem^^'^^ 
that the first-order change in energy due to a small per- 
turbation of the Hamiltonian is equal to the change in 



band structure energy, keeping the potential and den- 
sity fixed. Different Hamiltonian matrices based on DFT 
can therefore be expected to contain all information pro- 
vided that they are constructed from the same potential 
and density. 

Using the matrices and c^, coupling elements can 
be obtained through the transformation 



Vak — / .(^a'gJ^a'k' ^k'k^ 
a'k' 

_ ^A* qAM M 
^ak — 2_^^a'a^a'k' ^k'k- 

a'k' 



(10) 

(11) 



Thereby we have all the parameters in the Newns- 
Anderson Hamiltonian (Eq. (|5])), although the adsorbate 
state has an overlap Sak — {ci\k) with each of the metal- 
lic eigenstates, meaning the basis set is non-orthogonal. 
Below we will use expressions derived by Grimley^^ for 
the non-orthogonal case. 

If the adsorption induces a change 5p{e) in the metallic 
density of states, the adsorption energy can be written 
as 



ads 



J —c 



Sp{e)ede - ANep + na{eF - (12) 



where ei? is the Fermi level. The first term is the con- 
tribution to the binding energy from hybridization with 
the adsorbate (the factor of 2 assumes that each spin hy- 
bridizes equally and independently). The integral of the 
induced density of states 



AN ■ 



Sp{e) de 



(13) 



is the change in number of electrons below the Fermi 
level. If this integral is nonzero, either too much or too 
little charge will be counted in the integration up to the 
Fermi level, and the extra or missing electrons must then 
be deposited onto or taken from the Fermi level. This 
electron count correction is the second term, ANep- Fi- 
nally if the atom contributes Ua electrons which come 
from the adsorbate level Ca, this amount of extra charge 
must in turn be deposited on the Fermi level cf (last 
term). In the Newns- Anderson model, the first two terms 
of Eq. (p!2|) are expressed as the integral over a function 
r]{e) such that 



where 



tan7^(e) 



r]{e) de^UaicF - ^a), 



A(6,6) 



A(6) 



A(e)' 



A(e,e' 



/de', 



A(e, e') = TT ^ \esak - Vak\'^S{e' - Ck)- 



(14) 

(15) 

(16) 
(17) 
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FIG. 15: (Color online) (a) A(e), A(e) and e - Ca for H on 
58-atom Au cluster. The dotted line indicates the Fermi level, 
(b) Projected DOS on atom in eV~^ and total DOS of isolated 
cluster in arbitrary units, (c) Cumulative induced DOS. 



The one-variable function A(e) = A(e,e) is referred to 
as the chemisorption function and plays the role of a 
continuous couphng matrix element. r]{e) is the phase 
shift of the complex self-energy A(e) — zA(e) and is related 
to the the induced density of states by 



1 dr]{e) 
TT de 



5p{e). 



The PDOS on the adsorbate can be written as 



1 



A(e) 



TT (e 



A(e))2 + A2(e)- 



(18) 



(19) 



Because of the approximations used in this method, 
calculated binding energies are by themselves not useful 
(or accurate) compared to the DFT results. The strength 
of this method lies in the conceptual simplification that 
the binding energy can be understood from continuous 
functions such as A(e) and r]{e). This allows the origin 
of the coupling and binding energy to be attributed to 
particular states in the cluster. Next we will apply this 
to H, O, Li and F on a 58-atom Au cluster to understand 
the effect of magic numbers on chemisorption. 



A. Adsorption of H 

We perform a DFT calculation on a 58-atom Au cluster 
with H adsorbed to obtain the Hamiltonian and overlap 



FIG. 16: (Color online) (a) A(e), A(e) and e - Ca for O po. or 
Py on 58-atom Au cluster. Ca is indicated on the energy axis, 
(b) Projected DOS on all O 2p states in eV~^ (solid line), 
on px and py only (dashed line), and total DOS of isolated 
cluster in arbitrary units, (c) Cumulative induced DOS. 



matrix. In this calculation, only the Is basis function is 
included on H, but otherwise the parameters are identical 
to those used in earlier calculations. For the Is basis func- 
tion we calculate A(e) and A(e) which is shown on Figure 
[T5k together with the line e — e^. When e — = A(e) 
and A(e) is small, there will be resonances in Pa{^) (as 
per Eq. (p!9|) ) corresponding to states on the atom. The 
adsorbate level (circle) and Fermi level (dotted line) 
are indicated. The resulting PDOS is shown on Figure 
[T5b together with the total DOS of the Au cluster. This 
reveals that it is the strong coupling to low-lying metallic 
states (e ~ — lleV) which gives rise to a bonding, local- 
ized state at the bottom of the Au s-band, at —12 eV, and 
an antibonding state consisting of several peaks mostly 
above the Fermi level. 

On a side note, the very high PDOS at the antibonding 
state may seem surprising. Given that the adsorbate level 
Ca ^ —9 is much closer to the bonding state, the bonding 
state would be expected to be similar to |a), and thus 
have a high PDOS on the atom, while the antibonding 
state should be more similar to the metallic states and 
therefore have a low PDOS on the atom. However the 
inclusion of overlap Sak in the model causes part of the 
states on the neighboring metal atoms to be counted in 
the adsorbate PDOS, contributing to the prominence of 
the antibonding peak. While the overlap affects the cal- 
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culated PDOS, the overlap is correctly taken into account 
in binding energies and other parts of the formalism. 

In the creation of bonding and antibonding states 
the original adsorbate state has been eliminated, and 
a change 5p{e) in DOS has been induced in the clus- 
ter. The cumulative induced DOS —r]{e)/7T is shown on 
Figure [T5b . While the newly created bonding state at 
— 12eV can accept a certain amount of charge, a similar 
amount of charge has been removed from the remainder 
of the cluster DOS (mostly around — lOeV) such that the 
total integral of the induced DOS up to the Fermi level 
is zero. The extra electron from the H atom is there- 
fore deposited on the Fermi level. A higher-lying Fermi 
level implies a weaker adsorption energy, since the elec- 
tron is deposited at a higher energy. This is why clus- 
ters just past a magic number, characterized by a higher 
Fermi level but an otherwise similar spectrum, adsorb H 
more weakly than clusters just before a magic number. 
We can also see how the induced DOS integrates to zero 
only because the Fermi level is located at a gap between 
electronic shells: Within each electronic subshell there 
are fluctuations in the induced DOS which correspond 
to slight movements of the electronic shells but without 
the introduction of any extra charge. These cause 7^(e) to 
locally deviate from 0. The adsorption energy therefore 
may not depend simply on the Fermi energy in general, 
but must do so at the magic numbers. These results are 
consistent with previous findings for very small clusters, 
that H atoms effectively contribute their electron to the 
LUMO, behaving like an extra Au atom^^'— 



2.5 I I T- n 1 1.2 




Energy [eV] 



FIG. 17: (Color online) Projected density of states (top) on 
adsorbate (eV~^) and cluster (arbitrary units), and cumula- 
tive induced density of states (bottom) for Li (left) and F 
(right). 



causing a change in binding energy opposite that for H 
adsorption as seen from the DFT calculations. 



B. Adsorption of O 

Since the Newns- Anderson model only takes a single 
state into account while O has three p-states, we will 
assume that each of the states hybridizes independently 
and contributes to the adsorption energy as per Eq. (p!4|) . 
Thus we consider one Hamiltonian of the form (jSj for 
each p-state with varying Vak and Sak- 

The 2px and 2p2y-states are close to degenerate and 
have almost identical chemisorption functions. Figure 
[T6k shows the average A(e) and A(e) from the O 2px 
and 2p^-states. In this case the weaker splitting leads to 
greater smearing of the states close to the d-band and 
between the electronic shells. The higher-lying peaks in 
A(e) correspond to coupling with of the electronic shells. 
Fi gure fTBb shows the total PDOS pxyz{^) due to all three 
p-states (full lines) along with the contribution pxy{^) 
from px and Py (dashed). The most profound feature is 
the state between the top of the d-band and the Fermi 
level, which therefore is filled. The induced DOS inte- 
grates to 3.0 at the Fermi level (Figure [T6b), allowing 
space for six electrons counting spin-degeneracy. Since 
only four electrons are contributed, a total of two elec- 
trons are taken from the Fermi level into available lower- 
lying states. An upward shift in Fermi level therefore 
implies that more energy is gained from this transfer. 



C. Adsorption of Li and F 

Li is the simplest of the four cases. Here the bonding is 
weak enough that no significant splitting occurs. The ad- 
sorbate state instead broadens into a resonance far above 
the Fermi level, see the left part of Figure [TJl without in- 
ducing any states below the Fermi level. The electron 
contributed by the Li atom therefore moves down to the 
Fermi level causing the same dependence of adsorption 
energy on Fermi level as for H. 

F couples more weakly than O, and the bonding states 
are therefore split up less than for O (right part of Fig- 
ure [T7j). Since both bonding and antibonding states are 
occupied, F behaves like O except only one electron can 
be transferred from the Fermi level, meaning that the 
change in adsorption energy at magic-number clusters is 
generally smaller than for O. 



D. Comparison to Ft clusters 

Finally we shall briefly consider the binding of O to the 
58-atom Pt cluster. Again the p^^ and Py states are close 
to degenerate, and the average of their chemisorption 
functions is shown on Figure [T8k . Two primary features 
appear in the chemisorption function: a strong coupling 
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FIG. 18: (Color online) Chemisorption of O on 58-atom Pt 
cluster in comparison with Au from Figure [TBI 

within the d-band around e = — 8eV, and a number of 
higher-lying peaks corresponding to electronic shells like 
those of Au clusters. Due to the broader and higher- lying 
d-band, the adsorbate state splits into peaks over a wider 
energy range as seen on FigurefTSb. The increase in bind- 
ing on Pt compared to Au manifests itself as an increase 
in area below the curve in Figure [T6b cf . Eq. (p!4|) . The 
overall upward shift of the coupling leads to an upward 
shift of the induced density of states, and so only approx- 
imately 2.2 out of the 3 states contributed by O are in 
this case located below the Fermi level. The partial oc- 
cupation of O 2p-states has been studied and confirmed 
experimentally.— Recall that for Au, essentially all of the 
adsorbate-induced states were located below the Fermi 
level (Figure [T6b). 

The s-electron shell structure is mostly visible in the 
chemisorption function well above the Fermi level. Fur- 
ther, since the location of the Fermi level within the d- 
band prevents abrupt changes in the Fermi level with 
cluster size, the s-electron shell structure — as expected 



— cannot exert a strong influence on the chemical bind- 
ing on Pt clusters. 



VII. CONCLUSIONS 

The structure of very small Au clusters is intricately 
dependent on the s-electron hybridization, with clusters 
at magic numbers having very large band gaps. Clusters 
with different numbers of atoms deform considerably to 
minimize the band structure energy by creating gaps at 
the Fermi level. For odd-numbered clusters this results 
in a singly-fllled state in the middle of the gap, causing 
strong even-odd oscillations of HOMO and LUMO. 

Clusters based on regular geometries or a simple EMT 
potential show a more clear electronic shell structure and 
have large band gaps only at the major shell closings. 
These structures are less realistic, but are computation- 
ally feasible to optimize for larger cluster sizes. 

Adsorption energies of atoms on regular Au clusters 
oscillate with the electronic magic numbers. While lo- 
cal geometry is known to be important, the variation in 
binding energy of O due to magic numbers alone may be 
up to 1 eV. Clusters just before or after magic numbers 
are found to exhibit roughly halogen-like and alkali-like 
behavior while magic-number clusters are, as expected, 
universally unreactive. 

A more detailed analysis attributes the increase or de- 
crease in binding energy of speciflc adsorbates at magic 
numbers to properties of the adsorbate-induced density 
of states. Adsorption of O or F induces states below the 
Fermi level, allowing the transfer of electrons from the 
Fermi level into the lower- lying states. In contrast H and 
Li, despite having very different adsorbate levels and elec- 
tronegativity, only induce states above the Fermi level, 
and the electron contributed by these atoms is therefore 
transferred to the Fermi level. 
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